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■ NASA goals for increased aircraft 
efficiency 

■ Proof of concept architecture for flight 
control with load feedback 

■ Optimal control allocation with load 
constraints 

■ Wind energy challenges 

■ Wind turbine blade damage 

■ Adaptive contingency control 



SUGAR Ray design by NASA sponsored team 
led by The Boeing Company 
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Subsonic Fixed Wing Project 




■ Reduce environmental impact of aviation 

■ Increase aircraft efficiency 

■ Improve mobility of aircraft in airspace 

■ Unacceptable community noise and other 
environmental emissions 

■ Need to reduce fossil fuel consumption 

■ Demands from NextGen airspace 

■ Air transportation plays key role in our 
economy and quality of life 

■ Create prediction and analysis tools for design 

■ Develop concepts and technologies for significant 
improvements in noise, emissions and performance 

■ Partner with academia, industry, and government 
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Weight Reduction for Increased Efficiency 



■ Create new fabrication processes for lightweight 
materials, esp. large structures 

■ Design lightweight wing structures with aeroelastic 
tailoring to eliminate heavy control surfaces 

■ Use aerogels for super-lightweight insulation 

■ Increase temperature capability of composites for 
greater use in engines 



Electron Beam Freeform 
Fabrication (EBF3) 
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Polymer-enforced aerogel 



Fan containment system with 
high temperature capability 








Formation flight for drag reduction 



D-8 (double bubble) has natural laminar 
flow on wing bottom, design by MIT team 


> Increase laminar boundary layer by delaying transition to turbulent flow 

■ Design surfaces with favorable pressure gradients (natural laminar 
flow) 

■ Include active or passive local suction surfaces (hybrid laminar flow) 

> Use active aeroelastic tailoring of wing to reduce drag during cruise 

> Advanced CONOPS - formation flight 

> Develop & validate CFD codes for design & analysis of advanced drag 
reduction concepts 
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NextGen N+3 Concept Aircraft 
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Traditional Flight Control System 



> Traditional Approaches for Control Allocation 

■ Ganging of Actuators: use elevator for pitch, ailerons for roll, rudder for yaw 

■ Mixers: fixed combination of surfaces to achieve commands 

> Structural Limits 

■ Design engineers determine critical load paths in aircraft 

■ Mostly concerned with bending, torsion, and shear loads 

■ Load limits are determined through ground tests and flight tests 

■ Load limits imposed by restricting flight envelope; position & rate limiting actuators 
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New Challenges for Flight Control Systems 



■ Many redundant effectors 

■ Surfaces affecting multiple axes 

■ Actuator rate & position saturation 

■ Low control authority 

■ Lighter more flexible structures 



Optimal Control Allocation 

Given B, a desired vector u pnd e>0,, find u such that 
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Feasibility Study 



Objective: 

Use multiple control surfaces in most effective way, while remaining 
within structural load limits 

Approach: 

Replace traditional control allocation with optimal control allocation 
with load constraints and real-time load feedback 

■ Measure internal (structural) loads along critical load paths 

■ Use aircraft aerodynamic and structural models to determine 
incremental loads due to incremental surface deflection 

■ Include structural load constraints and measured loads in optimal 
control allocation problem 

Significance: 

This approach enables fuel efficient aircraft with many multi-purpose 
control surfaces to achieve acceptable performance & safety 
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Study Assumptions 


■ Only considering static loads due to lift and rolling forces 

■ Finite element model (FEM) of aircraft wings and tail 

■ Loads due to lift and roll are applied to nodes in FEM model 

■ Bending moments calculated using finite element analysis (FEA) 

■ A select number of load points are monitored and included in the 
optimal control allocation constraints 
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Proposed Framework 



a d 

desired accelerations 

U P 

preferred delta surface positions 

B 

control effectiveness matrix 

T 

incremental loads matrix, where Tu gives the 
incremental loads at critical points 

y s 

structural loads from sensors (or model in sim) 

F E 

external forces due to lift and body moments 

F, 

internal structural loads at critical points 
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Optimal Control Allocation 



Optimal Control Allocation Problem: 

Given B, a desired vector u d and £>0, find u such that 


J - Bu— a 


d 


+ 8 


u — u 


d 


V 

error 

min 




V 

control 

min 


00 

J 


is minimized subject to: 

u m ^<u<u 
where F I max are critical point load limits 
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Can also be formulated with load minimization, using 
constraints given above and cost function: 
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Up-scale GTM Simulation 



5.5% GTM-T2 Solid Model 


Up-scale GTM Simulation 

■ Simulink model based on 5.5% dynamically 
scaled aircraft derived from wind tunnel & flight 
test data 

■ Up-scaled by incorporating Reynolds adjusted 
aero tables 

■ Actuator models sized for up-scale GTM 

■ NASA Glenn’s Simp2 engine (simplified version 
of C-MAPSS40k) 

■ GTM bare airframe 

■ 6 ailerons, 4 elevons, 2 rudders, 2 stabs, 2 flaps 

■ Vehicle Management System 

- sensor processing module 

- mission manager 

- guidance/control 

■ Vehicle Control Augmentation System 

- reference model dynamic inverse controller 

- optimal control allocator 
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Finite Element Modeling Approach 



■ Wings/tails modeled as cantilever beams with fixed ends at roots 

■ Constant thickness hollow aluminum shells following outer mold 
line give beam cross section properties 

■ Beam nodes located at centroids of wing cross sections 



Vortex-lattice model of aircraft Wing cross sections 
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Finite Element Model of GTM 


■ Beam mesh for each wing has 20 nodes and 1 9 beams 

■ Each beam has 6 degrees of freedom - 3 translation & 3 rotation 







U. Wyoming 9/27/11 


FEM of GTM aircraft 


FEM of wings and empennage 




Finite Element Analysis 



■ Assume static loads and static 
response, F = Kx 

■ Stiffness matrix K is derived 
from FEM 

■ K _1 is computed off-line 

■ Static loads applied during 
simulation to FEA nodes 

■ Measured loads are calculated 
from deflections using K A 
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Application of External Static Loads 



■ Elliptically distributed lift load applied to 
nodes along wings and horizontal tail 

■ Loads arising from roll moments applied as 
concentrated forces in z-direction on each 
aileron in proportion to aileron deflection 
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Integration with Simulation 


■ Flap-wise bending moment at critical points are 
calculated & passed to control allocator 

■ K' 1 and B are used during simulation to determine 
incremental loads matrix T 

■ Aileron forces are assumed to be proportional to 
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Roll forces 
applied here 
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Flight conditions: 

■ Altitude 30,000 ft 

■ Mach 0.85 

Test cases: 

■ Normal case 

Load limits set to values 
determined for safe flight 

■ Case I 

Right aileron load limit set to 
55,000 ft- lb 

■ Case II 

Left outboard aileron deflection 
limits set to ±0.01 deg 
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Surface deflection (d) 


Test Cases 


Example of control 
allocation gone bad 

COMMANDED ROLL MANEUVER 

Both produce same roll moment 

Structural load limits exceeded 

f M I M ^ - 

Structural loads within limits 
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Yaw Rate (d/s) Pitch Rate (d/s) Roll Rate (d/s) 


Roll Doublet Case I 


Normal Load Constraints 




o.i 1 1 — — — i 1 

-0 1 1 1 1 1 1 

0 2 4 6 8 10 
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Bending Moment (ft-lb) Bending Moment (f-lb) 


Bending Moment for Right Outboard Aileron 



Case I: Outboard aileron critical point limit set to 55,000 ft-lb 



01 23456789 10 


Normal 

Load 

Constraints 


Reduced 

Load 

Constraints 


Time (s) 
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Right Aileron 3 Right Aileron 2 Right Aileron 1 

Deflections (d) Deflections (d) Deflections (d) 



Aileron Deflections Case 1 





Reduced load limits 


U. Wyoming 9/27/2011 


22 




Right Aileron 3 Right Aileron 2 Right Aileron 1 

Deflections (d) Deflections Deflections (d) 


Aileron Deflections Case li 




Restricted position limit 
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Formulation & Solution Approaches 


Mixed Optimization Formulation 


Find u that minimizes J = II (CB)u-a d 
subject to u <u<u , |w|<zi 

J min max’ max 


+ £ 


u-u p 


❖ 






Solution Approaches Using Different Norms 

norm: ||Hli=^KI 


Convert to linear programming problem 

- Simplex algorithms (Bodson) 

- Interior-point algorithms (Peterson, Bodson) 


4 norm: 




- Active Set Method with norms squared (Harkegard) 

- Interior-point algorithms 

4, norm: IHL^maxCw,-) 


- Simplex algorithm (Bodson, Frost) 
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Right Aileron Deflections (d) Left Aileron Deflections (d) 


Normal Case with L versus /, Norm on Control 


L Norm on Control 


4 Norm on Control 
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Right Aileron Deflections (d) Left Aileron Deflections (d) 


Normal Case with L versus l x Control Norm 



L Control Norm /, Control Norm 
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Conclusions & Future Work 


Proposed framework performed adequately in simulation & 
proof of concept demonstration was successful 

■ Control allocation 

- Try load minimization 

- Use weights on surfaces and critical points depending on health of components 

- Explore non-feasible solutions 

■ Loads model 

- Investigate robustness, computation time, sensitivity 

- Include torsion, structural dynamics 

■ Work with sensors to measure loads 

■ Include aeroelasticity in simulation 

■ Flight test technology at Dryden Flight Research Center 
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Our Next Challenge: Distributed Control Effectors 



Objective: 

• Develop variable stiffness materials for distributed control effector skins 

Approach: 

• Investigate mechanisms that can impart variable stiffness to material 
systems to enable novel control effectors to control lightweight flexible 
wings 

• Test various concepts on distributed control effector model to determine 
feasibility as a structural element 

Significance: 

• Variable stiffness material systems can enable control of aeroelastically 
tailored lightweight wings to meet SFW fuel burn goals 

Results: 

• Bench model (shown top right) developed to study angles and 
deformations of distributed control effectors to develop requirements 
and test candidate variable flexibility control surface skins 

• This activity will involve materials, structures, aeroelasticians, controls 
and dynamics experts working concurrently to design, analyze, build 
and test a distributed control surface concept 




Above: Bench model of distributed control surface 
Below: Sketch of aeroelastically tailored wing with 
distributed control effector on trailing edge 
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Wind Industry Observations 


Wind Industry Practices 

> Industry is relatively low-tech and 
very protective of IP 

> Research funding is limited 

> Very little vertical stratification 

> New technologies need quick 
and cost-effective integration 
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Wind Industry Challenges 

> Building large turbines (>5 MW) 

> Operating & maintenance costs 

> Turbine reliability 

> Grid integration 

> Community noise 

> Wind farm siting 



Condition Monitoring in Wind Turbines 



■ SCADA system 

> Supervisory Control and Data Acquisition for wind farm 

> Medium- and long-term changes in environmental & 
operating conditions 

> Minimal fault diagnosis 

> Lots of data, not always useful 

■ Short-term condition monitoring 

> Equipment set up for one month for vibration, acoustic, 
strain, nacelle acceleration testing 

■ Scheduled maintenance & inspection 

■ Acceptance of CM by operators/developers 

> Dependent on cost of CM system 

> Might affect warranty 

■ Some OEMs are moving towards guaranteed uptime 



Image: www.vertigo.net.au 
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Leading Causes of Blade Failures 1 



Shear webs 
(glass over core) 


1) Manufacturing defects - wrinkles in laminate, missing or incomplete bond lines, dry fibers 

2) Progressive damage initiating from leading-edge erosion, skin cracks, transport, 
handling, or lightning strikes 

3) Excessive loads from turbine system dynamics or dynamic interaction with control system 

4) Out-of-plane forces and distortion of blade sections (“bulging/breathing” effect) mostly in 
root transition region, due to blade loading 

5) Excessive loads due to unusually severe atmospheric conditions 

^NV Renewables, Seattle, WA, “Lessons Learned from Recent Blade Failures: Primary Causes and Risk-Reducing 
Technologies”, D.A. Griffin & M.C. Malkin, 49 th AIAA Aerospace Sciences Meeting, Jan 2011, paper no. 2011-259 
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Flexible Structure Control 



Flexible Structure Control Challenge: 

■ Structural modes can be excited by feedback control 

■ Low pass filter can reduce problems, but they have limitations 

■ Residual mode filter (RMF) has model of structural mode, including 
phase and frequency, that can be removed from feedback 

■ Flexible structures ARE intrinsically modal systems 


■ Flexible aerospace structures are difficult to model and they operate 
in poorly known environments 

■ Adaptive control helps, but requires minimum phase plants (ASPR) 

■ Residual Mode Filters (RMF) can cancel transmission zeros, restoring 
ASPR 


Recall: 


A system( A, B, C ) is ASPR whenCB > 0 and its closed- loop 
t ran sf erf unction P( s ) = C(sl — A) 'Bis minimum phase. 
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Partition Plant into ASPR & non-ASPR 


Assume original system (A ,B ,C ) can be partitioned into: 
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Adaptive Controller using RMF 
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Nonlinear Wind Turbine Plant 
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Adaptive Pitch Control for FAST Simulator* 


■ Objective: 

■ Input: 

■ Output: 

■ Disturbance: 

■ How: 

■ RMF 


Generator speed for 
turbulent wind input 

— Baseline PI 

— Adaptive RMF 


Regulate generator speed and reject disturbances 
Rotor speed 

Collective blade pitch, constant generator torque 
Step function 

Model gusts as step functions 

Designed for drive-train rotational flexibility mode 



U. Wyoming 9/27/2011 


*NREL’s FAST simulator of CART2 (high fidelity simulation of flexible 2-bladed wind turbine) 
see: http://wind.nrel.gov/designcodes/ 
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FAST Turbine Blades 


FAST blade configuration files: 

■ 21 distributed stations along span 

■ Flapwise & edgewise stiffness 

■ Flapwise & edgewise bending modes 

Assumption: 

Blade damage can be 



Department of Wind Energy, 
Risoe National Laboratory 


Flapwise 


Edgewise 
'rection 


represented by reduction in flapwise & edgewise stiffness 


Damaged blade configuration files: 

■ Edgewise and flapwise stiffness are varied at a blade station 

■ Blade bending mode shapes are recomputed 

■ Structural damping and other parameters were left unchanged 
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Blade Node Sensitivity to Stiffness Changes 



Full factorial study performed to determine node sensitivity: 

■ Parameters: blade damage, wind speed, blade pitch 

■ Levels: 8 for damage, 7 for wind, 1 0 for blade pitch 


Loads on blades 
are primarily due 
to aerodynamic 
forces 



X deflection vs Damage Level vs Node 


Damage(%) 


Node 


1 0.325 

1 0.32 

1 0.315 
Ko.31 


- 0.305 

- 0.3 

- 0.295 

- 0.29 


0.285 

0.28 


■ 0.275 
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Study of Turbine Response to Blade Damage 


Preliminary study of effects of blade stiffness reduction 

■ Damage located on one blade at station 7, 30% from blade root 

■ Study run in open-loop with no generator a 

speed tracking 

■ Generator torque held fixed at rated 
torque 

■ Simulation run with steady wind speeds 
from 12-24 mps 

■ Collective pitch varied from 0. 1 -0.45 
radians 

■ Blade tip displacement was measured 
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Tip Displacement vs Pitch vs Damage 
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Collective Pitch Residual Stiffness Fraction 


■ Pitch & wind speed dominate change in stiffness 

■ Damage detection tool needs to factor out impact of pitch and 
wind speed to use deflection 
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Standard Deviation of 
Out-of-plane Tip Deflection (m) 


Change in Tip Deflection with Power Output 



Hypothesis: Reducing power output through generator set-point 
reduction will reduce loads on turbine blades 



Out-of-plane tip deflection std. dev. for 3 generator set-points & 7 damage levels. 
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De-rating Generator for Reduced Loads 



Time (s) 

Simulation Wind Input 
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Top: Generator set-point 
Bottom: Generator speed 



Control Integrated with Condition Monitoring 



> Simulation demonstrating contingency controller lowering generator 
set-point for turbine with blade damage when winds are turbulent & 
above rated speed 

> Resulting decrease in blade root bending could extend service life 


Blade root bending moment 
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